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Summary 

A pulsed thermocouple is used for measuring gas 
temperatures above the melting point of common 
thermocouples. This is done by allowing the 
thermocouple to heat until it approaches its melting 
point and then turning on the protective cooling gas. 
This method requires a computer to extrapolate the 
thermocouple data to the higher gas temperatures. In 
earlier work by this author the extrapolation was 
done by using a first-order exponential curve fit to 
predict the final thermocouple wire temperature. 
Since radiation effects were neglected, the gas 
temperature was not computed. Hand calculations 
had to be used to estimate the gas temperature. This 
report describes a method that includes the effect of 
radiation in the extrapolation. Computations of gas 
temperature are provided, along with the estimate of 
the final thermocouple wire temperature. Results 
from tests on high-temperature combustor research 
rigs are presented. 


Introduction 

An earlier investigation by the author (ref. 1) 
described the use of a pulsed thermocouple to 
measure gas temperatures above the melting point of 
common thermocouples. This method of measuring 
temperature is intended for the measurement of 
temperatures at the exit of experimental aircraft 
combustors at temperatures to 2400 K and pressures 
to 4 MPa (40 atm). The previous investigation 
described an approach that uses a thermocouple 
cooled by a small jet of inert gas. When a 
measurement is to be made, the cooling jet is turned 
off and the thermocouple is allowed to heat up to 
near its melting point. When the temperature of the 
thermocouple approaches its melting point, the 
cooling is reapplied. The data are then fitted to a 
first-order exponential function. The final 
temperature that the thermocouple would have 
attained is then calculated by extrapolation. 

The computer program (ref. 1) did not take into 
account the fact that at the higher temperatures the 
heating curve deviates from a true exponential. This 
deviation is the result of radiant energy (obeying 
Stephan’s T^ law) being absorbed or emitted by the 
thermocouple wire. 

The analysis described in this report takes into 
account the T* radiation terms in the differentitil 


equation describing the temperature of the 
thermocouple wire as a function of time. The report 
describes the solution of this differential equation for 
time as a function of temperature. This solution 
cannot be inverted (except numerically) to give 
temperature as a function of time. A computer 
program is described that fits measured data to the 
theoretical curve based on this more complete 
analysis. The computer program uses the gradient- 
expansion method (ref. 2) to fit the data to the 
theoretical function. The program computes final 
thermocouple wire temperature and final gas 
temperature. 

This report also presents typical input and results 
for the computer program. Data and results are 
discussed from tests in two combustor test facilities. 


Theory 

This section describes the theoretical equations 
necessary to compute gas temperatures with a pulsed 
thermocouple. Most of the time the thermocouple is 
protected with a jet of cooling gas, as shown in figure 
1. When a temperature measurement is to be made, 
the cooling gas is turned off and the thermocouple 
output is sampled at a high rate and recorded. Just 
before the thermocouple reaches its melting point the 
cooling is reapplied to protect the thermocouple wire. 
The gas temperature can then be calculated by 
extrapolation from the initial heating curve. For the 
extrapolation to be valid, it must be based on a 
theoretical heating curve. The derivation of the 
theoretical equation is described here. All symbols 
are defined in appendix A. 

The equation that describes the pulsed- 
thermocouple wire temperature can be derived from 
the basic heat transfer relations (ref. 3). Assume a 
bare wire thermocouple with infinitely long leads in a 
hot gas stream. This assumption causes the 
conduction effects to be neglected. Very little error is 
introduced if we neglect the transfer of heat to the 
junction by conduction along the wire for carefully 
designed probes. Thus in the absence of conduction, 
heat can be transferred to the wire by convection of 
the gas, by radiation from the gas, and by radiation 
from the duct walls. Also heat can be transferred 
away from the wire by radiation. 

The rate of heat storage in the wire will be equal to 
the rate of heat entering the wire minus the rate of the 
heat leaving the wire. The rate of heat storage Qs per 
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practice the Psc must be determined experimentally 
and generally falls in the range 0.8 to 1.0. 

The rate of heat transfer by radiation is given by 
(ref. 3) 


?, = <7e,,[(l-a)7j + eg7^-7tt]7rZ» (4) 

where a is the Stefan-Boltzmann constant, is the 
emissivity of the wire, a is the effective absorptivity 
of the gas, is the duct temperature, and eg is the 
emissivity of the gas. The first term in equation (4) 
represents the heat received by the wire from the hot 
walls of the duct. The second term represents radiant 
heat received from the gas. The third term represents 
radiant heat emitted from the wire. 

Combining equations (1) to (4) gives 


dt^ 


-KidT^ 

Tt+K2T^-K^ 


(5) 


Figure 1. - Thermocouple heating curve. 


where 


unit length is given by 


Qs ~Qc 


( 1 ) 


^1 = 


pCD 

Aae^ 


( 6 ) 


where is the rate of heat convected per unit length 
to the wire by the gas and ^^is the net heat radiated 
per unit length to the wire. 

The rate of heat storage per unit length of the wire 
is given by (ref. 3) 


Ki = 


'NuKgPsc 

Daey, 


and 


(7) 


Qs=pC 


ttD^ dTy/ 
4 dt 


(2) /T3=/:2 7g + [(l-«)7^-Feg7^' 


( 8 ) 


where p is the wire density, C is the specific heat of 
the wire, is the wire temperature, t is the time, 
and D is the wire diameter. 

The rate of heat transfer to the wire by convection 
qc is given by (ref. 3) 


To solve equation (5), we integrate both sides of 
the equation. The integration is easier if we factor the 
denominator. The roots of a fourth -order equation 
can be found by algebraic methods (ref. 4). The roots 
of the equation are 


= 7T Nu Kg Psc {Tg - Tyy) (3) 

where Nu is the Nusselt number. Kg is the thermal 
conductivity of the gas, Psc is the probe shape 
constant, and Tg is the gas temperature. For an 
infinitely long wire in crossflow Psc is unity. The 
probe shape constant was introduced to take into 
account the fact that the presence of a probe to 
support the wire will cause a reduction in the 
effective Nusselt number of the thermocouple. In 


Tyy=ai ±iP,a2,a^ (9) 

where 
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Equation (5) can then be rewritten as 


( 12 ) 

(13) 


(14) 


E= —2/32 (2a 1 _q ;2 _Q.j) 
and 

E’=2[/3(ai -a2)(ai -a3)-/32] 
Thus 7^3 can be rewritten as 


7/3 = -^+/ 


£■ 2 +/^ £• 24 ./;^ 

=^3A +1H3B 

Equation (16) can then be integrated to get 
r = //lln(a2 - 7’w) + 7/2ln(rH,-a3) 


( 21 ) 


( 22 ) 


(23) 


dt = 


-Ki dT^ 

(JTy,-a\ -i/3)(7’w-«l +'^)(7’w -«2)(^w -«3) 

(15) 


+ //3Aln[(7'^-«l)2+)S2] 

+ 2//3Btan - • ( ) + //4 (24) 


or 


dt = 


Hx H 2 

7^-02 Ty,-oii 








(16) 


where 


Hx = 


-Kx 

(a2 -a3)[(a2 -«l)^ +/3^1 


(17) 


H 2 = 


-Kx 

(as -a2)[(o!3 -«l)^+ /32] 


(18) 


H3 = 


-Kx 

(ai - tt 2 + i/3)(ax - a 3 + //3)(2/d) 


(19) 


If the denominator of Hs is multiplied out, we get 


Hs = 


-Kx 

E+iF 


( 20 ) 


where 


where H 4 is a constant of integration. 

Equation (24) shows the theoretical relationship 
between the wire temperature Tyy and the time i. In 
general all the parameters in the equation are known 
except for the gas temperature Tg, the probe shape 
constant P^c, and the integration constant H 4 . After 
a measurement a set of wire temperature readings are 
known. The procedure used finds the values of Tg, 
^SCf and H 4 that result in the best fit of the 
temperature data to the theoretical equation (eq. 
24)). The next section describes the computer 
program written to fit equation (24) to the data. 


Description of Computer Program 

The FORTRAN IV computer program described in 
this report is designed to calculate gas temperature by 
using data taken from a separate pulsed- 
thermocouple controller. A listing of the program 
and its various subroutines is shown in appendix B. 
The program input requirement is a set of wire 
temperatures taken at regular time intervals, the 
Mach number, the total pressure, the wall 
temperature, and the probe shape constant. The 
computer program output is the extrapolated wire 
temperature and the computed gas temperature. In 
addition, if the probe shape constant has not been 
entered, the computer program will calculate and 
output PSC, the probe shape constant. 

The program uses a curve-fitting procedure from 
reference 2 called the gradient-expansion method to 
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fit the theory to the input data. Two parameters, gas 
temperature tgas and possibly probe shape constant 
PSC, are Adjusted for best fit of the theory to the 
data. These parameters are adjusted until the sum of 
the squares of the differences between the measured 
wire temperature and the theoretical wire 
temperature is a minimum. The error, which is called 
CHiSQR, is defined by 


n 2 

CHISQR= ^ |^(^data)y “ (^theory)jj (25) 

where Tdata is the measured wire temperature, 
^’theory is the corresponding theoretical wire 
temperature, and n is the number of measured data 
points. Note that the theoretical wire temperatures 
must be evaluated point by point at the same values 
of the time parameter used for the measured data. 

Both the gradient-expansion procedure and the 
evaluation of CHiSQR require computation of 
theoretical wire temperature at every measurement 
time. In addition, the gradient-expansion method 
requires values for dT^ldTg and bTy/lbP^c at every 
measurement time. These requirements create a 
difficulty because the analytical solution to the 
differential equation expresses time as a function of 
wire temperature in equation (24). The equation 
cannot easily be inverted to yield the needed wire 
temperature as a function of time and its derivatives. 
As a result a great amount of the computer time is 
devoted to numerically inverting the equation and 
evaluating the derivatives. Since theoretical wire 
temperature values at the measurement times are not 
available directly from equation (24), they are 
calculated by interpolating in a table of wire 
temperature-time pairs that do satisfy equation (24). 
This table must be regenerated whenever equation 
parameters are changed. 

This procedure must be repeated once for every 
evaluation of wire temperature and twice for every 
evaluation of the derivatives. The derivatives are 
approximated by computing the differences in wire 
temperature that result for two values of the 
parameters TGAS and PSC: one value slightly above 
the present value and one value slightly below the 
present value. 

The main computer program takes care of reading 
the input data, calling the curve-fitting routines, 
deciding when the curve fit is good enough, and 
writing the results. Initially input data of Mach 
number, pressure, and duct temperature are read in 
as well as 1000 readings of thermocouple wire 
temperature. The temperatures represented by these 
numbers are taken at equal time intervals before and 
during the temperature rise. The first 100 readings 


represent the thermocouple wire temperature while 
the cooling air is on. The rest of the 900 temperature 
readings are taken during the temperature rise of the 
thermocouple wire when the cooling air is turned off. 
If the cooling air is turned on again before the 9(X) 
readings are taken, the remaining readings are zero. 

After the data are read in, a call to subroutine 
STCFIT determines the best estimate of the 
temperature ramp starting time. This is necessary 
because the theoretical curve is always forced to pass 
through this point. 

With the starting time determined, the curve- 
ntting process begins. Repetitve calls to curfit and 
FDERiv result in adjustments to several parameters 
such that CHISQR is decreased. With every 
adjustment in the parameter values a call to congen 
is needed to evaluate the constants in equation (24). 
The parameters adjusted include the gas temperature 
tgas; the probe shape constant PSc; and flamda, a 
parameter whose value controls the curve-fitting 
process. The probe shape constant is adjusted only if 
its value is not included in the input data. If the PSC is 
to be adjusted, the variable nterms is set equal to 2 
by the computer program; otherwise NTERMS is set 
equal to 1 and only tgas is adjusted. Thus the main 
program recalls fderiv and curfit until the decrease 
in CHISQR is less then 1 percent. This value of 1 
percent was chosen by trial-and-error methods to 
provide a wire temperature within 1 or 2 K of the 
ultimate wire temperature without using an 
unreasonable amount of computer time. 


Subroutine curfit 

Subroutine CURFIT makes a least-squares fit to a 
nonlinear function by using the gradient-expansion 
algorithm described in appendix C. The algorithm is 
really two curve-fitting techniques combined into one 
program. One of the techniques works well when the 
variables are far from the correct values, and the 
other works well when they are close to the final 
values. A parameter X (called flamda in the 
program) is used to change the curve-fitting routine 
gradually from one technique to the other. 

The subroutine works by starting with 
FLAMDA = 0.001 (when flamda is less than 1 the 
fitting technique that works close to the minimum is 
dominant — see appendix C). The error (appendix 
C) between the measured and theoretical data is 
called both CHiSQl and CHISQR in the program. 
CHISQI is an initial value of calculated once when 
the subroutine is entered. The program makes 
changes in the wire temperature, the probe shape 
constant, and flamda until a new value of x^ (called 
CHISQR) starts to decrease, at which time flamda is 
divided by 10 and the subroutine returns to the 
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calling program. It is the responsibility of the calling 
program to check CHISQR to see if the change in 
CHiSQR since the last call to CURFIT is small enough 
to stop the program. If it is not, subroutine CURFIT 
should be called again without changing the value of 
the current flamda. 


Subroutine FDERiv 

Subroutine fderiv computes data needed by the 
curve-fitting routine curfit. The data needed are the 
derivatives of the wire temperature with respect to 
both gas temperature and the probe shape constant. 
Also needed are theoretical values of wire 
temperature evaluated at the measured time (the 
times corresponding to the measured wire 
temperatures). The derivatives are determined from 
(ref. 5) 

dT^ r^(Tg -F 1 , PSC) - T„(Tg - 1 , PSC) 

dTg 2 ^ ^ 

dTy, _ Ty,(Tg. PSC -F 0.001) -7’w(7’g. PSC-0.001) 
Opsc “ 2(0.001) 

(27) 


If the probe shape constant is not to be calculated 
(NTERMS=l)> only equation (26) will be calculated. 
The theoretical values of wire temperature are 
generated from equation (24) with a call to 
subroutines tabl and intrp. 

The subroutine returns a 1000- by 3-element eirray. 
The derivative of the wire temperature with respect to 
the gas temperature at time / is returned in array 
DERiV(i,i). The derivative of the wire temperature 
with respect to the probe shape constant at time / is 
returned in array deriv(I,2). The table of the 
computed wire temperatures at time / is returned in 
array DERiV(i,3). 

Function xicalc 

Function XICALC computes the sum of the squares 
of the differences between the measured wire 
temperature and the theoretical wire temperature 
(from the numerically inverted equation (24)). The 
sum of the squares of the differences will be 


RANGE 2 

XICALC = Yt [rw(/)-7’theory(7)l (28) 

1= START I J 


The program first calls subroutine CONGEN to 
generate new constants for equation (24) since the gas 
temperature and the probe shape constant may have 
changed. Subroutine tabl is then called to generate a 
table of theoretical temperatures and times. The 
interpolation necessary is done by this subroutine 
and not by subroutine INTRP because the output of 
this routine is a single number, the error xicalc, and 
not an entire table of numbers. 

Subroutine tabl 

The purpose of subroutine tabl is to generate 
values of theoretical wire temperatures and times for 
subroutine intrp. Subroutine CURFIT, fderiv, and 
function XICALC require a value of theoretical wire 
temperature at every measurement time. These wire 
temperatures must be obtained by inverting equation 
(24). However, because of the form of equation (24) 
a numerical inversion will have to be done. A call to 
subroutine tabl generates a table of temperature- 
time pairs that satisfy equation (24). Then a call to 
INTRP interpolates in this table to get temperatures at 
the measurement times. 

To generate the interpolation table, a set of 
temperatures is needed to put into equation (24) to 
obtain computed times. The values of computed time 
that result from equation (24) should be as close as 
possible to the measured times for accurate 
interpolation by subroutine INTRP. The set of 
temperatures is determined one at a time, starting 
with a known point on the theoretical curve. Each 
succeeding temperature is computed from the 
previous one by using a linear approximation to the 
theoretical curve (fig. 2). The linear approximation 
will have a slope equal to the slope of the theoretical 



Figure 2. - Graphical representation of linear approximation to 
theoretical wire heating curve. The tj are measured times 
and the tj are computed times from equation (24). 
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curve at the previous temperature. Thus each 
succeeding temperature will be 


'y + i 


= 7’. -I ^ 


Tj 


(29) 


where j= 1,2,3,..., n measured data points. The times 
corresponding to the measured data points are tj. 
The times tj are computed by evaluating equation 
(24) with = Tj. The derivative of equation (24) is 


dt _ -Hi H2 

dTw ai-T^ T’w-aa 



Time 


, 2[//3A(T>,-ai )-//3B/3] 

(7’,-«i)2+^2 ( > 

In the program tj+\ -t'j is defined as DELTIM and 


DELTMP = 


DELTIM 

{dt/dTy,)j 

* ‘J 


(31) 


The program starts by setting Tj = T\ = tave, which 
is the temperature on the theoretical curve; and 
t'j=t{ is equal to mstime-start. The next 
temperature + i is evaluated by setting 
O + l =^2 =MSTIME*(start+ 1) in equation (29). 
What results is a table of theoretical time- 
temperature pairs that do satisfy equation (24), 
where the times are not exactly equal to the 
measurement times. The array of times is called 
TIMC, and the array of temperatures is called TC in 
the program. A linear interpolation will need to be 
done because temperatures at the exact measurement 
times are needed. 


Subroutine intrp 


Figure 3. - Search process for subroutine STCFIT. 


the cooling is turned off. This line is called tave. The 
other line is the best fit through approximately the 
first 50 points of the temperature rise. Since a 
solenoid is used to turn the cooling air on and off, 
there will be some delay between when the power is 
removed and when the cooling air actually stops 
flowing. The solenoid power is turned off at data 
point 100, and the starting point search ranges 
between data points 100 and 130. 

The starting point of the search process is shown in 
figure 3. A standard least-squares fit to a straight line 
of the data from point 100 to point 160 is performed. 
In general, point 100 is not the true starting point; so 
this line Qine 1 in fig. 3) will not intersect the tave 
line at point 1(X). In fact, if the starting point of the 
data for the least-squares line is varied from 100 to 
130, the intersection of the least -squares line Qine 2) 
with TAVE will approach the true starting point and 
then back away. Therefore the intersection point will 
have a maximum as the starting point is varied. The 
output of this routine is this maximum value of the 
starting point. This represents the best 
approximation to the start of the ramp. 


Subroutine INTRP is used to correct the table of 
theoretical temperatures (array TC) generated by 
subroutine TABL. Subroutine intrp performs a linear 
interpolation between the calculated data points so 
that the calculated times (and corresponding 
temperatures) fall exactly on the measured time. The 
resulting interpolated values of temperature are 
stored in array TC. 


Subroutine CONGEN 

Subroutine CONGEN computes the constants 
necessary to evaluate equation (24). Constants K \ , 
K 2 , and K-j are evaluated by using equations (6) to 
(8). The wire emissivity for clean platinum was 
found to be (ref. 6) 


Subroutine stcfit 

Subroutine stcfit determines the starting point of 
the thermocouple temperature rise. The starting 
point is defined as the intersection of two straight 
lines. One line is the best fit through the data before 


e„,»0.085 + (0.76E-4)r^/ (32) 


where Tyff is the final wire temperature in K. The 
other parameters used for platinum (type R) 
thermocouple wire are (ref. 7): 
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Wire density, kg/m^ 0.2078 x 10^ 

Stefan-Boltzmann constant, 

J/K'* sec m2 0.56697 X 10~2 

Wire specific heat, J/kg K 0. 1427 x 10 

Gas effective absorptivity 0 

Gas effective emissivity 0 

Wire diameter, m 0.8128x 10 


Kg = (0.3007 E-3 ).tgasO -^8 j/(sec K m) (33) 

Nu = 188.41 ♦ (VwDI A» MN • P) • TG AS “0 ® 

• [l +0.2*(MN)2| 

where wdia is the wire diameter, mn is the Mach 
number, p is the pressure in pascals, and is in K. 
The gas effective absorptivity and emissivity are 
assumed to be zero. This corresponds to a 
transparent gas and the worst case for radiation 
effects. 

The subroutine also computes aj , «2, 03, /3, , 

Hi, /^3A. and //3B from equations (10) to (23). The 
value of is computed by putting the initial 
conditions into equation (24) and solving for 
The initial temperature is the average cooled 
temperature tave. The initial time is the 
measurement time interval mstime times start. 

Function evaltm 

Function evaltm evaluates equation (24) to 
obtain a calculated time for an input of wire 
temperature. The input wire temperature must be 
between the initial average cooled temperature tave 
and «2 in order to avoid taking the logarithm of a 
negative number. Values of oq, ai, 03, /3, H\, Hi, 
Ht,p^, and H^ must have been previously calculated 
with a call to the congen subroutine. 

Subroutine matinv 

Subroutine matinv does an inversion of a 1- or 
2-degree matrix. For a 1 -degree matrix only a simple 
reciprocal is needed. For a 2-degree matrix the 
adjoint matrix is calculated. Then each element is 
divided by the determinant to form the inverse 
matrix. The original matrix is then replaced by its 
inverse. 


Tests and Results 

A pulsed-thermocouple system was tested in a 
combustor rig at the Ah Force Wright Aeronautical 
Laboratory (AFWAL) as part of a joint af-nasa 


program on instrumentation. The system included a 
probe, a sample-and-hold voltmeter, a 
microcomputer-based controller, and a digital 
recorder, as shown in figure 4. Figure 5 shows the 
probe that was put into the combustor. The probe 
consisted of a water-cooled shell with a replaceable 
platinum (type R) thermocouple. Compressed-air 
cooling for the thermocouple was controlled by a 
fast -acting solenoid valve. The thermocouple voltage 
was converted to digital form by a sample-and-hold 
digital voltmeter. A microcomputer was used to 
control the voltmeter and turn the cooling air on and 
off. The time between data points (called mstime) 
was controlled at 0.0042 second. This value was 
chosen so that most of the ramp would be included in 
the 1000 data points. If a different probe with a 
different time constant were used, this mstime would 
have to be changed. 

A full curve including the final wire temperature 
could be recorded for each pulse because the gas 
stream of the combustor configuration under test 
was not hot enough to require the cooling air to come 
on. The data were first processed by the computer 
program to compute the probe shape constant. The 
average computed probe shape constant for 20 pulses 
at fixed combustor conditions was 0.91, with a 
maximum deviation of 0.09. This deviation is the 
result of the fact that the burning process is not 
constant during the pulse and thus results in a 
temperature that can vary during the pulse by as 
much as 3.2 percent. 

With the average probe shape constant of 0.91 the 
data were curve fit 60 percent of the way up the 



Figure 4. - Block diagram of pulsed-thermcouple system. 
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Figure 5. - Thermocouple probe. 



Time, sec (scale times 0. (XM20) 

Figure 6. - Resuits with final temperature well below wire 
melting point. 


curve. It is estimated that at least 60 percent of the 
curve could be measured at the highest expected gas 
temperatures. A typical result is shown in figure 6. 
The solid line is the measured data (a total of 1000 
data points). The triangles and squares represent the 
theoretical curve. The squares represent the portion 
of the curve that was used in the computation. The 
triangles represent the portion of the curve that was 
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Figure 7. - Results with final temperature near wire melting 
point 


extrapolated by using the theoretical curve. The 
computed final wire temperature varied from 1525 K 
to 1581 K, with an average of 1561 K for the 20 
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readings. The actual final wire temperature varied 
from 1525 K to 1575 K because of fluctuations in the 
burning. A comparison between the final wire 
temperature computed using 60 percent of the ramp 
and the actual final wire temperature measured for 
the 20 readings showed a maximum deviation of 3 
percent. 

The average of the 20 computed gas temperatures 
was 1691 K, with a maximum deviation of 47 K, or 
2.7 percent. The difference of 130 K between the 
computed wire temperature and the gas temperature 
is the radiation error. It is estimated that the 
radiation error can be computed to within about 20 
percent, which for this case would be ±26 K. 

Results for a pulsed-thermocouple probe different 
from the probe just described were obtained during a 
high-temperature combustor test at the Lewis 
Research Center as shown in figure 7. The probe 
shape constant for this geometry was determined at 
lower temperatures than shown in figure 7 to be 0.96. 
The gas temperature for the data shown in figure 7 
was 2300 K, and the final computed wire temperature 
was 2190 K. The wire melts at 2215 K. The protective 
compressed air was set to turn on at about 2000 K in 
order to assure a long thermocouple life. 


Concluding Remarks 

The pulsed thermocouple was developed as an 
instrument to determine high gas temperatures. The 
pulsed feature is needed at temperatures above the 
melting point of common thermocouples or when 
streaking of a combustion process is occurring. The 
cooling gas was found to adequately protect the 
thermocouple during this high-temperature 
operation. 

The computer program for computing gas 
temperature was designed to take the radiation 
error into account. The program requires as input the 
Mach number, the wall temperature, and the total 
pressure in addition to the thermocouple data. Tests 
at temperatures below the melting point of platinum 
thermocouples show that the pulsed-thermocouple 
system can compute the gas temperature to within 
about 4 percent with as little as 60 percent of the 
temperature step as input data. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, December 15, 1980 
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Appendix A 
Symbols 

Definition 

symbol 

symbol 


a 


parameter of function 

C 

SPHT 

specific heat of wire 

D 

WDIA 

wire diameter, m 


H1,H2,H3A, 

intermediate constants 


H3B.H4 


Kb 


thermal conductivity of gas, J/sec K m 

Kx,K^K^ 

K1.K2.K3 

intermediate constants 


MN 

Mach number 

Nu 

NU 

Nusselt number 

P 

P 

pressure. Pa 

P 

^SC 

PSC 

probe shape constant 

Qc 


rate of heat transferred by convection into surface of wire, J/sec m 

Qr 


rate of heat transferred by radiation, J/sec m 

Qs 


rate of heat storage, J/sec m 

T 


temperature, K 


TAVE 

average temperature 

Td 

TDUCT 

duct temperature 

Tg 

TGAS 

gas temperature 

Ty, 

TWIRE 

wire temperature 

Twf 

TWF 

final wire temperature 

t 


time 

X,Y 


general independent variables 

Yl 


intermediate constant 

a 

ALPHAG 

effective absorptivity of gas 

ai,a2. 

ALPHAl, 

intermediate constants 

«3 

ALPHA2 



ALPHA3 

BETA 

intermediate constant 


EGAS 

emissivity of gas 

f»v 

E1+E2*T 

emissivity of wire 

a 

SIGMA 

Stefan-Boltzmann constant, J/K^ sec m^ 

x2 

CHISQR, 

least-squares error 


CHISQl 
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Appendix B 
Computer Programs 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c; 

c 

c 

c 


c 


c 

c 

c 

c 


c 


c 


c 


ROUTINE FOR CURVE FITTING DATA FROM A PULSED THERMOCOUPLE, 
THERMOCOUPLE DATA SHOULD BE CONSTANT FOR THE FIRST 100 
DATA POINTS (THERMOCOUPLE COOLED), THE THERMOCOUPLE IS THEN 
HEATED ON AN EXPONENTIAL HEATING CURVE (900 DATA POINTS), 
THE PROGRAM NEGLECTS CONDUCTION ERRORS, 


INPUTS REQUIRED AREJ 

WIRE TEMPERATURE (1000 DATA POINTS) (K), 

MACH NUMBER, 

PRESSURE (Pa), 

D U C T T E M P E RAT U RE ( K ) , 

PROBE SHAPE CONSTANT, 

ALSO THE FOLLOWING PARAMETERS MUST BE SET TO 
THE PROPER VALUE DEPENDING ON THE TYPE OF 
T H E R M 0 C 0 LI P I... E U S E D J 

MSTIME := TIME BETWEEN MEASUREMENTS (THIS PROGRAM), 
WDIA WIRE DIAMETER, (SUBROUTINE CONGEN), 

WDENS == WIRE DENSITY, (SUBROUTINE CONGEN), 

SPHT - WIRE SPECIFIC HEAT, (SUBROUTINE CONGEN), 
WIRE EMISSIVITY, (SUBROUTINE CONGEN), 

EGAS = EMISSIVITY OF GAS, (SUBROUTINE CONGEN), 
ALPHAG = ABSORPTIVITY OF GAS, (SUBROUTINE CONGEN), 


REAL TW I RE ( 1 ()()() > f TA VE . ALPHAl r ALPHA2 v ALPHAS f BETA ? 

1 HI r H2 , H3A r H3B r H 4 » TDUCT i- MSTIME r MN r P ? MNN . MTMP 

INTEGER START » RANGE 

REAL TC ( 1 000 , 2 ) r T I MC ( 1 000 f2)f DER I V ( 1 000 3 ) 

DATA MSTIME/0, 42E-2/ 

MSTIME IS IN SECONDS, 

INTEGER IrNTERMS 

REAL CH I SQO f FL AMDA TGAS r TWF f PSC r CH I SQR r X 

COMMON /BLK 1 /TW IRE » TA VE . ALPHAl ALPHA2 r ALPHAS y BETA y 
1 HI y H2 y H3A y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 

COMMON /BLK2/ST ART .RANGE 
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COMMON /BLK3/TCrTIMCrDERIV 
C 

c 

c 

C READ TEMPERATURE DATA (TWIRE) 

C IN DEG, K , (1000 DATA POINTS), 

C 

DO 20 1=1 r 1000 
20 READdrSO) TWIRE(I) 

READ INPUT DATA 

WRITE (7. 30) 

30 F0RMAT( IX f ' INPUT MACH NUMBER') 

READ ( 5 r 40) MN 
40 F0RMAT(F7,4) 

WRITE(7r50) 

50 FORMAT( IX. ' input PRESSURE IN Pa,') 

READ(5r60) P 
60 F0RMAT(F11 ,3) 

WRITE(7? 70) 

70 FORMAT( IX r ' INPUT DUCT TEMPERATURE IN DEG, K,') 

READ(5f80) TDUCT 
80 F0RMAT(F9,2) 

WRITE(7f 90) 

90 FORMAT < IX r ' INPUT PROBE SHAPE CONSTANT,') 

READ (5 >100) PSC 

100 F0RMAT(F6,3) 

THE NEXT 3 STATEMENTS ARE NEEDED ONLY 
FOR THE EXAMPLE IN THIS REPORT, 

WRITE(7> 101 ) 

101 FORMAT ( IX > 'MACH NUMBER TEMPERATURE DEG, K,') 
READ(5>80) MTMP 

AVERAGE COOLED WIRE TEMPERATURE, 

TAVE = 0,0 
DO 110 1=1 >99 
10 TAVE = TAVE+TWIRE(I) 

TAVE = TAVE/99, 

DETERMINE START OF RAMP, 

115 CALL STCFIT 

C 

C TEMPERATURE OVER MELTING POINT? 

C 

118 DO 120 RANGE=100> 1000 

IF (TWIRE(RANGE) ,LE,400, ) GO TO 130 
120 CONTINUE 

12 
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130 RANGE = RANGE-1 

CUROE FIT. 

CHISQO ••= 0. 

FLAMDA == 0.001 
TWF = TUI RE (RANGE) 

TGAS = TUF 
NTERMS =-- 1 

1 35 I F ( PSC . NE . 0 . ) GO TO 140 
NTERMS = 2 
PSC 0.9 

■FDERIO" COMPUTES THE DERIOITIOE OF TWIRE UITH RESPECT 
TO TGAS & PSC. ALSO IT RETURNS OALUES OF CALCULATEfi 
THEORETICAL WIRE TEMPERATURE AS A FUNCTION OF TIME. 

140 CALL FDERIV(TGAS!-PSCnNTERMSrTUF) 

■CURE IT- MODIFIES TGAS AND PSC TO OBTAIN THE BEST- 
MATCH BETWEEN THE THEORETICAL CUROE AND THE ACTUAL DATA. 

CALL C U R F I T ( N T E R M S » P S C r T G A S » C H I S Q R y F 1.. A M D A r T W F' ) 

•CHISQR* IS THE ERROR BETWEEN THE THEORETICAL CURUE AND 
THE ACTUAL MEASURED DATA. IF THERE IS LESS THAN A 
ONE PERCENT CHANGE IN THE ERROR SINCE THE LAST 
CALL TO CURFIT THEN THE PROGRAM IS FINISHED. 

X = ABS( (CHISQR-CHISaO)/CHISQR) 

IF (X.LT.0.01) GO Ta 150 
CHISQO = CHISQR 
TWF = ALPHA2 
GO TO 140 

150 WRITE(7yl60) TGAS 

160 FORMAT ( IX f "GAS TEMPERATURE = 'yF9.2y' K") 

WRITE(7rl70) ALPHA2 

170 FORMAT( IXr "FINAL WIRE TEMPERATURE - "yF9.2y" K " ) 
WRITE(7ylS0> PSC 

180 FORMAT ( IX y "PROBE SHAPE CONSTANT = "yF6.3) 

STOP 123 
END 


SUBROUTINE CURFIT ( NTERMS y PSC y CHISQR y FLAMDA y TWF ) 
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c purpose: 

C THIS SUBROUTINE MAKES A LEAST SQUARES CURVE FIT TO 

C A NON-LINEAR FUNCTION. 

C 

C TIME = SET OF INTEGERS TAKEN AS INDEPENDENT VARIABLE. 

C TWIRE = ARRAY OF WIRE TEMPERATURE READINGS TAKEN 

C AS DEPENDENT VARIABLE. 

C START = INTEGER VALUE OF TIME FOR START OF DATA. 

C RANGE INTEGER VALUE OF TIME FOR END OF DATA. 

C NTERMS = NUMBER OF PARAMETERS (MAX. == 2). 

C TGAS == PARAMETER It GAS TEMPERATURE. 

C PSC = PARAMETER 2t PROBE SHAPE CONSTANT. 

A =: ARRAY OF PARAMETERS. 

FLAMDA - PROPORTION OF GRADIENT SEARCH INCLUDED. 

TWF = ESTIMATED FINAL WIRE TEMPERATURE. 

CHISQR CHI SQUARE FOR FIT. 

SUBROUTINE CURE IT (NTERMS ^ PSC f TGAS » CHISQR f FLAMDA r TWF ) 


REAL TWIRE ( 1.000 .1 r TAVE i- ALPHA! y ALPHA2 y ALPHAS y BETA y 
1 H! y H2 y HSA y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 

INTEGER START y RANGE 

REAL TC( 1000y2) y TI MC ( 1000 y 2 ) y DERI V ( ! 000 y 3 .> 

REAL BE(2.1 yAL(2y2) y PSC y TGAS y CHISQ 1 y CHI SQR y FLAMDA y TWF y 
1 A(2) yB(2) yARRAY(2y2.) 

INTEGER I y J y K y NTERMS y ERFLAG 


COMMON /BLKl/TWIRE y TAVE y ALPHA! y ALPHA2 y ALPHAS y BETA y 
! H ! y H2 y H3A y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 

COMMON /BLK2/STARTy RANGE 

COMMON /BLK3/TCyTIMCyDERIV 


AL = ALPHA MATRIX. 

BE = BETA MATRIX. 

10 DO 20 J=!y NTERMS 

BE(J) = 0. 

DO 20 K=!yJ 
20 AL(JyK) := 0. 

TRUNCATE TGAS SINCE SMALL CHANGES IN TGAS 
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C CAUSE UNNECESSARY ITERATION 

C 

Ad) = AINT(TGAS) 

A(2) == PSC 

30 DO 60 I=START»RANGE 

DO 50 J=1»NTERMS 

BE(J) = BE(J)+(TWIRE(I)-DERI0(I>3))*DERI0(If J) 
DO AO K=1»J 

40 AL(JrK) = AL ( J i- K ) +DERIO ( I » J ) *DERI0 ( I f K ) 

50 CONTINUE 

60 CONTINUE 

DO 70 J:=;L i-NTERMS 
DO 70 K=lrJ 

70 AL(KfJ) = AL (J.K) 

C 

C EOALUATE CHI SQUARE AT STARTING POINT 

C 

CHISQl == XICALC(Ad ) >A(2) fERFLAGr TWF) 

C 

80 DO 100 J==1tNTERMS 

DO 90 K==l rN TERMS 

CALCULATE ALPHA PRIME MATRIX (CALLED ARRAY) 

AND ALSO INOERT IT» 

ARRAY( Jj K)==AL( JrK) 

90 ARRAY ( J r J ) ==ARRA Y ( J ? J ) * (1 ♦ +FLAMDA ) 

100 CONTINUE 

C A I- 1... M A T I N 0 ( A R R A Y y N T E R M S ) 

110 B(2) - A(2) 

DO 130 J--=- l rNTERMS 
B ( J ) = A ( .J ) 

DO 120 K==-l vNTERMS 

120 B(J) “ B(J) T BE(K)*ARRAY( J.K) 

130 CONTINUE 


TRUNCATE B(l) & B(2) TO CONSIDER ONLY INTEGER OALUES 
OF TEMPERATURE AND ONLY 2 SIGNIFICANT FIGURES FOR PSC. 

B(l) == AINT(Bd)) 

B(2) - B(2)*100. 

B(2) -- AINT(B(2)) 

B(2) = B(2)/100. 

CALCULATE CHISQR FOR NEW PARAMETER OALUES. 

CHISQR = XICALC(B(1 ) .B(2) .ERFLAG.TWF) 

ERFLAG=6 IF ALPHA2 IS TOO LOW. 

IF (ERFLAG.EQ.6) GO TO 140 
IF (CHISQl -CHISQR) 140.150.150 
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C 

C IF CHISQF^ INCREASED? INCREASE FLAMDA* 

C 

140 FLAMDA = lO.OKFLAMDA 

GO TO 80 
C 

C IF CHISQR DECREASED? DECREASE FLAMDA & SET NEW 

C VALUES FOR TGAS S PSC* 

C 

150 TGAS = B(l) 

PSC = B<2) 

FLAMDA = FLAMDA/10, 

160 RETURN 

END 


C 

C S U B R 0 U TINE F D E R I V ( T G A S ? P S C ? N T E R M S ? T W F ) 

C 

C 

C PURPOSE 

C COMPUTE THE DERIVATIVE OF TWIRE WITH RESPECT TO BOTH 

C TGAS & PSC AND ALSO EVALUATE THE THEORETICAL EQUATION. 

C 
C 

C TGAS = ESTIMATED GAS TEMPERATURE (K). 

C PSC PROBE SHAPE CONSTANT. 

C NTERMS = NUMBER OF TERMS (1 OR 2). 

C TWF = ESTIMATED FINAL WIRE TEMPERATTURE (K). 

C 

C COMMENTS 

C THIS PROGRAM COMPUTES THE DERIVATIVE OF TWIRE WITH 

C RESPECT TO TGAS AND STORES THE VALUES IN THE ARRAY 

C DERIV(I?1) WHERE 1=1? 1000. THE DERIVATIVE OF TWIRE 

C WITH RESPECT TO PSC IS STORED IN DERIV(I?2). THE 

C THEORETICAL EQUATION EVALUATED AT EACH MEASURED TIME 

C I:MSTIME*(TIME integer)! is stored in DERIV(I?3). 

C 

C 

C 

SUBROUTINE FDER I V ( TGAS ? PSC ? NTERMS ? TWF ) 


REAL TWIRE ( 1000) ?TAVE?ALPHA1 ?ALPHA2? ALPHAS? BETA? 
1 H 1 ? H2 ? H3A ? H3B ? H4 ? TDUCT ? MST I ME ? MN ? P ? MNN ? MTMP 
C 

INTEGER START? RANGE 
C 


16 



nono no ooo noon noon nno 


REAL TC<1000r2) »TIMC(J.000»2) ,DERH^( 1000y3) 

C 

REAL DELTA(2) » T , PC » TGAS » PSC r X » Y » Z » TUF 
C 

INTEGER FrFFfERFLAGrlrLl 
C 

DATA DELTA/.! tOj 0.001/ 

C 

COMMON /BLKl/TWIRE » TAOE r ALPHA! » ALPHA2 r ALPHAS » BETA r 
! H ! r H2 » H3A f H3B r H4 r TDUCT i- MSTIME » MN i- P » MNN f MTMP 
C 

COMMON /BLKlVSTARTrRANGE 
C 

COMMON /BLK3/TCTTIMCrDERIM 
C 
C 

C COMPUTE DATA FOR DERIOATIOE OF TWIRE WITH 

C RESPECT TO TGAS. 

C 

C COMPUTE DATA POINTS FOR T=TGAS+ DELTA ( ! ) 

C AND PC=-PSC. 

C 

!0 T = TGAS + DELTA(l) 

PC = PSC 

GENERATE CONSTANTS. 

CALL CONGEN ( T y PC ? ERFL AG y TWF ) 

IF (ERFLAG.NE.O) STOP 997 

GENERATE A TABLE OF COMPUTED TIMES AND TEMPERATURES. 

TABL PUTS DATA INTO TC AND TIMC 

CALL TABL<!) 

COMPUTE DATA POINTS FOR T=TGAS-DELTA ( ! .'i 
AND PC=PSC. 

20 T =•■ TGAS - DELTAd) 

PC = PSC 

GENERATE CONSTANTS. 

CALL CONGEN ( T y PC y ERFLAG y TWF ) 

IF (ERFLAG.NE.O) STOP 996 

GENERATE ANOTHER TABLE. 

CALL TABL (2) 

INTERPOLATE BETWEEN DATA POINTS SO THAT THE 
CALCULATED TIMES (TIMC) CORRESPOND TO THE 
MEASURED TIMES (I*MSTIME). 
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C 

C 

c 

c 

c 

c 


30 


35 


C 

C 

C 


INTERF'OLATE FOFi; T=:TGAS+DELTA (1. ) ♦ 
CALL INTRPdfl) 

INTERPOLATE FOR T=TGAS-DELTA ( 1 ) , 


CALL INTRP<2i-.1.) 


CALCULATE DERIVATIUE OF TWIRE WITH RESPECT TO 
TGAS AND STORE IT IN DERIO(Iyl)* 

DERIU ( START » .1. ) •--- 0. 

LI START + I 
DO 35 I=LIyRANGE 

DERI V ( I y 1 ) :::: ( TC < I y .1. ) --TC ( I y 2 ) ) / ( 2 . *DELTA CD) 
CONTINUE 

IF (NTERMS,EQ» .1. ) GO TO 60 

COMPUTE DATA FOR DERIMATIOES OF TWIRE WITH 
RESPECT TO PSC* 

COMPUTE DATA POINTS FOR T^TGAS AND 
PC=PSC+DELTA(2) ♦ 


T =: TGAS 

PC = PSC + DELTA(2) 


GENERATE CONSTANTS. 

CALL CONGEN ( T y PC y ERFLAG y TWF ) 

IF (ERFLAG.NE.O) STOP 995 

GENERATE A TABLE OF COMPUTED TIMES AND TEMPERATURES. 
CALL TABLd) 

COMPUTE DATA POINTS FOR T=TGAS AND 
PC=PSC-DELTA(2) . 


T = TGAS 

PC = PSC - DELTA (2) 

GENERATE CONSTANTS. 

CALL CONGEN ( T y PC y ERFLAG y TWF ) 
IF (ERFLAG.NE.O) STOP 994 

GENERATE ANOTHER TABLE. 

CALL TABL(2) 
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c 

C INTERPOLATE BETWEEN DATA POINTS FOR PC=PSC+DELTA ( 2 ) * 

C 

CALL INTRP(.ty2) 

C 

C INTERPOLATE FOR PC=PSC-DELTA ( 2 ) 

C 

CALL INTRP(2y2) 

C 

C CALCULATE DERIVATIVE OF TWIRE WITH RESPECT 

C TO PSC AND STORE IN DERIV(If2). 

C 

50 DERIV(STARTr2) 0. 

DO 55 I L .1.y RANGE 

DER I V ( 1 . 2 ) := ( TC ( I r ;t ) -TC ( I r 2 ) > / ( 2 * *DELTA ( 2 ) ) 

55 CONTINUE 

C 

C GENERATE A TABLE OF ONLY THE FUNCTION (TWIRE VS* TIME). 

C 

C GENERATE CONSTANTS FOR TGAS 8 PSC. 

C 

60 CALL CONGEN ( TGAS y PSC. ERFLAGfTWF) 

IF (ERFLAG.NE.O) STOP 993 
C 

C COMPUTE A TABLE. 

C 

CALL TABLCL) 

C 

C INTERPOLATE BETWEEN DATA POINTS 

C 

CALL INTRP( .1. .3) 

C 

C STORE THE INTERPOLATED FUNCTION (TWIRE VS. TIME) 

C INTO DERIV(Iy3). 

C 

70 DO 75 I:=.1. V 1.000 

DERIV(Iy3) := TC(I.l) 

75 CONTINUE 

RETURN 
END 


C 

C 

C 


C 

C 


FUNCTION XI CALC ( TGAS y PSC y ERFLAG y TWF ) 


PURPOSE 

TO COMPUTE CHI SQUARE FOR PRESENT PARAMETER VALUES. 
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C 

C 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 


c 


TGAS = ESTIMATED GAS TEMPERATURE <K), 

PSC = PROBE SHAPE CONSTANT. 

TWF = ESTIMATED FINAL WIRE TEMPERATURE (K). 


REAL TWIRE( 1000) r TAUE » ALPHAl ? ALPHA2 j ALPHAS r BETA ? 
1 HI j H2 r H3A y H3B r H4 . TDUCT y MSTIME y MN y P y MNN y MTMP 

INTEGER START y RANGE 

REAL TC( 1000y2) y TIMC ( 1000 y 2 ) y DERI U ( 1 000 y 3 ) 

REAL XyYyTCMyER 

INTEGER ERFLAG y J y I y L 1 y L2 y K 


COMMON /BLK 1 /TWIRE y TAME y ALPHAl y ALPHA2 y ALPHAS y BETA y 
1 HI y H2 y H3A y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 

COMMON /BLK2/STARTy RANGE 

COMMON /BLKS/TCyTIMCy DERIM 


GENERATE CONSTANTS FOR THEORETICAL EQUATION. 

10 CALL CONGEN(TGASyPSCyERFLAGyTWF) 

IF (ERFLAG.NE.O) RETURN 

GENERATE TABLE OF THEORETICAL DATA POINTS OF 
TEMPERATURE MS. TIME. THE MALUES OF TIME ARE ONLY 
APPROXIMATELY EQUAL TO THE ACTUAL MEASURED TIMES. 

CALL TABL(l) 

J = START 
LI = START+1 
L2 = RANGE+1 
ER := 0, 

INTERPOLATE SO THAT TIMES FOR THEORETICAL DATA 
CORRESPOND TO DATA FOR ACTUAL MEASURED TIMES. 

55 DO 90 I^LlyRANGE 

X = MSTIME*FLOAT( I > 

DO 60 K=JyL2 

IF (X.LT.TIMC(Ky 1 ) ) GO TO 70 
60 CONTINUE 
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70 J = 

IF NO CHANGE IN TC(Jfl) NO INTERPOLATION NECESSARY. 

IF (TC(Kfl).EQ.TC(Jfl)) GO TO 85 
TCM = TC(Jfl) + (X-TINC<Jyl))5fc(TC(K?l)-TC<Jfl))/ 

1 (TIMC(Ki-l)-TIMC(Jf .1.) ) 

CALCULATE XI SQUARED. THIS IS INCLUDED INSIDE 
INTERPOLATION LOOP FOR CONUIENENCE. 

80 ER = ER+(TWIRE(I>-TCM)**2 

GO TO 90 

85 TCM = TC(Jrl) 

GO TO 80 
90 CONTINUE 

XICALC=ER 
RETURN 
END 




ii 


C 

C 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 


c 

c 


SUBROUTINE TABL(F) 

PURPOSE 

GENERATES A TABLE OF THEORETICAL TEMPERATURE VS. 
TIME DATA POINTS. THE 'COMPUTED TIME CORRESPONDS 
CLOSELY WITH THE MEASURED TIME BUT NOT EXACTLY. 

F = SECOND INDEX ON TC AND TIMC. (F 1. OR 2). 


REAL TWIRE( 1000) , TAVE r ALPHA 1 r ALPHA2 y ALPHAS r BETA r 
1 HI f H2 y H3A y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 

INTEGER START y RANGE 

REAL TC( 1000y2) y T IMC ( 1 000 y 2 ) y DERI V (1000 y 3 ) 

REAL DELTIM y DELTMP y TCALC y ZY 
INTEGER FyJyIyl...lyL2 

COMMON /BLK 1 /TWIRE y TAVE y ALPHA 1 y ALPHA2 y ALPHAS y BETA y 
1 HI y H2 y H3A y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 
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c 

c 

c 

c 

.1.00 


COMMON ./BLK2/STARTff«ANGE 
COMMON / B L K 3 / T C f T .1 M C r n E R ;i: V 

TC( START.. F) == TAOE 
TCALC = TAVE 

MST.TME == ACTUAL TIME BETWEEN MEASURED DATA POINTS 
OF TWIRE, 

DELTIM == MSTIME 


FIX FIRST POINT, 

TIMC(START.F) EOALTM ( TCALC ) 

INITIALIZE POINT COUNTER FOR ACTUAL MEASURED TIMES, 

J = START 
Z = ALPHA2-0,01 

COMPUTE TABLE OF TEMPERATURES 

.10 LI = START+1 
L2 == RANGE+1 
DO 140 I=:L1.L2 

COMPUTE TCALC AT POINT I - 1. 

IF TCALC > Z WE ARE ON TOP FLAT PORTION OF CUROE 
AND TEMPERATURE WILL NOT CHANGE ANY MORE, 

IF (TCALC, GT,Z) GO TO 130 

COMPUTE DERIOATIOE OF TIME WITH RESPECT TO 
TEMPERATURE FOR THEORETICAL CURVE, 

Y = -H1/(ALPHA2-TCALC) + H2/ ( TCALC-ALPHA3 ) + 

1 2 , 0* ( H3A* < TCALC-ALPHAl ) -H3B*BETA ) ./ 

1 ( (TCALC-ALPHAl )**2+BETA**2) 


DELTIM ■■= CHANGE IN TIME FROM THE LAST DATA NECESSARY 
TO MAKE THE CURRENT DATA POINT FALL APPROXIMATELY ON 
THE ACTUAL MEASURED TIME, 

DELTMP = THE CORRESPONDING CHANGE IN TEMPERATURE, 
NOTE THAT THE THEORETICAL FUNCTION (EVALTM.) GIVES 
TIME BACK FOR AN INPUT OF TEMPERATURE, 
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DELTMP = DELTIM/Y 
TCALC = TCALC+DELTMP 
TC(IrF) = TCALC 
TIMC(IyP) = EVALTH( TCALC) 

C 

C DETERMINE THE CHANGE IN TIME NECESSARY FOR THE NEXT 

C CALCULATED DATA POINT TO FALL ON THE NEXT MEASURED 

C DATA POINT* 

C 

120 J == J+1 

DELTIM = MSTIME3t!FL0AT( J+1 )-TIMC( I f F) 

IF (DELTIM.LE.O.O) GO TO 120 
GO TO 140 
C 
C 


130 

J = Jfl 



TC<IyF) =- 

ALPHA2 


TIMCd yF) 

- MSTIME*FLOAT( J) 

140 

CONTINUE 



RETURN 



END 







c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 


c 

c 

c 


c 


c 


c 

c 


SUBROUTINE INTRP( F y FF ) 


PURPOSE 

TO INTERPOLATE BETWEEN DATA POINTS CALCULATED FROM 
THE THEORETICAL EQUATION* USED BY SUBROUTINE FDERIU* 


F = SPECIFIES THE SECOND INDEX (1 OR 2) FOR TC AND TIME* 
FF ••== SPECIFIES THE SECOND INDEX (lr2r OR 3) FOR DERIU* 

UPON RETURN TC(I-F) HAS THE INTERPOLATED VALUES. 


SUBROUTINE INTRP<FyFF) 


REAL TWIRE( 1000) y TAVE y ALPHA 1 t ALPHA 2 y ALPHAS y BETA r 
1 HI y H2 y H3A y H3B y H4 y TDUCT . MSTIME y MN y P y MNN y MTMP 

INTEGER START y RANGE 

REAL TC( 1000y2) yTIMC< 1000 y 2) y DERI V ( 1 000 y 3 ) 
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REAL X 


C 

INTEGER F?FFj J-r I f KrLl. 

C 

COMMON /BLKl/TWIRE » TAOEi- ALPHA 1 1 - ALPHA2 r ALPHA3 . BETA . 
1 HI y H2 r H3A y H3B y H4 y TDUCT y MST I ME y MN y P y MNN y MTMP 
C 

COMMON /BLK2/STARTy RANGE 
C 

COMMON /BLK3/TCyTIMCyBERIV 


10 LI == STARTfl 


160 

170 

175 


180 


C 


1 

1 


C 


190 


200 


DERIM(IyFF) HAS NOT BEEN USED YET AND 
IS USED AS TEMPORARY STORAGE 

DO 170 I=--lylOOO 
DERIVd.FF) = TC(IyF) 

J = START 

DO 200 I-U...1 ..RANGE 

X --- MSTIME>l(FLOAT<I) 

K J~1 
K = KTl 

IF (X,GE.TIMC(K.F) .> GO TO 180 
J = K-1 

IF (DERIO(KyFF) .EQ.DERIU< JyFF) .> GO TO 190 

TC(IyF) = DERIM(J.FF) T ( X-T I MC ( J . F ) .) * 

( D E R I M ( K y F F ) ■ - D E R 1 0 ( J y F F ) ) ./ 
(TIMC(K.F.L-TIMC( JyF) .> 

GO TO 200 

TC(IyF) DERIO(JyFF) 

CONTINUE 

RETURN 

END 


C 

C 


C 

C 

C 

C 

C 

C 


SUBROUTINE STCFIT 


PURPOSE 

TO CALCULATE THE STARTING LOCATION OF THE RAMP, THIS 
LOCATION IS CALLED "START". 
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nn on on nan nnnn on n n n a nnon 


C INPUT IS TEMPERATURE-: DATA "TWIRE"* THE DATA POINTS 

C 1 = 100 TO 160 WII..L BE CURVE FIT TO A STRAIGHT LINE 

C Y = M*X -f B WHERE Y=TWIRE(I) AND X=I=TIME. 

C 

C 

SUBROUTINE STCFIT 

THIS ROUTINE CALCULATES THE STARTING LOCATION 
OF THE RAMP, 

REAL TWI RE ( 1 000 ) , TA VE » ALPHA 1 » ALPHA2 f ALPHA3 f BETA » 

1 HI j H2 » H3A » H3B f H4 » TDUCT r MSTIME » MN » P » MNN r MTMP 

INTEGER START f RANGE 

INTEGER IrJrK^IX 

REAL SUMX r SUMY y SUMXX y SUMYY y SUMXY y B y M y X 

COMMON /BLKl /TWI RE y TAUE y ALPHA 1 y ALPHA2 y ALPHA3 r BETA y 
1 HI y H2 y H3A y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 

COMMON /BLK2/STARTy RANGE 


10 START = 100 

USE K TO INCREMENT STARTING INDEX FOR TWIRE FROM 
100 TO 130, 

DO 50 K=100yl30 

INITIALir/E VARIABLES 

SUMX = 0, 

SUMY = 0, 

SUMXX = 0, 

SUMYY = 0, 

SUMXY = 0, 

J = 0 

PERFORM STANDARD LEAST SQUARES CURVE FIT TO A 
STRAIGHT LINE, 

20 DO 30 I=Kyl60 

J = J-fl 

SUMXY = SUMXY-f (FLOAT( I ) )*TWIRE( I ) 

SUMY = SUMY-fTWIRE( I ) 

SUMX = SUMX+FLOAT( I ) 

30 SUMXX = SUMXX-f (FLOAT(I) )»*2 

CURVE FIT DATA TO Y = M*X f B 
C 


25 



on nnnonnnrjnnonn noonnnnnnnn 


M= ( ( FLOAT < J ) ) *SUMXY--SUHX)kSUMY ) / ( ( FLOAT ( J > ) *SUMXX-SUMX*SUMX ) 
B (SUMY-M*SUMX)/<FLOAT( J> ) 

STRAIGHT LINE CURVE FIT DONE* 


DETERMINE IF LAST X WAS A MAXIMUM? 

IF NO CONTINUE. 

IF YES PROGRAM IS DONE AND 'START" IS SET EQUAL 
TO X AS THE STARTING LOCATION OF THE RAMP. 

SET Y = TAVE 

40 X (TAVE-B)/M 

IX = IFIX(X) 

IF ( IX. GT. START) START = IX 
50 CONTINUE 

RETURN 
END 


SUBROUTINE CONGEN < TGAS f PSC r ERFL AG ^ TWF ) 

PURPOSE 

THIS PROGRAM GENERATES THE CONSTANTS NEEDED FOR THE 
THEORETICAL TIME VS. TEMPERATURE EQUATION (EQU. 24). 

TGAS = ESTIMATED GAS TEMPERATURE (K). 

PSC - PROBE SHAPE CONSTANT. 

ERFLAG == AN ERROR FLAG => SET=0 IF NO ERROR. 

TWF == ESTIMATED FINAL WIRE TEMPERATURE. 


SUBROUTINE CONGEN ( TGAS r PSC r ERFLAG r TWF ) 


REAL TWIRE (1000 ) r TAVE f ALPHAl r ALPHA2 r ALPHA3 » BETA » 

1 HI rH2f H3ArH3BrH4f TDUCTrMSTIMErMNf pjMNNf MTMP 
C 

INTEGER START f RANGE 
C 

REAL NU r K1 r K2 f K3 r PSC » WDI A r WDENS , SPHT » S I GMA r 
1 El fE2f EGASf ALPHAGfXr AZ»ErBZfF»Y 
C 

DATA WDIA/0.8128E-3/» WDENS/0 . 20785E+5/ » SPHT/0 . 1 427E+3/ i- 
1 SIGMA/0. 56697E-7/r El/0.85E~l/» E2/0.76E-4/r 
1 EGAS/0. O/r ALPHAG/0.0/ 
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C 

c temperatuf?e: == deg* k* 

C WDIA = WIRE DIAMETER (METERS)* 

C WDENS = WIRE DENSITY ( K£Jtri/m**3 ) * 

C SF'HT = WIRE SPECIFIC HEAT ( J/ ( Ksim » K ) ) * 

C SIGMA = STEPHAN BOL.TZMAN CONSTANT ( J/ ( SEC * y K**4 » it.**2 ) ) * 

C EMISSIVITY OF WIRE = E1.+E2)KTWF* NO UNITS ON El* 

C E2 HAS UNITS OF 1/DEG* K* 

C EGAS = EMISSIUITY OF GAS* 

C ALPHAG = ABSORPTIVITY OF GAS* 

C P = PRESSURE (PASCAL)* 

C 

C 

INTEGER ERFLAG 

AZyBZyErF ARE TEMPORARY VARIABLES* 

AZ 8 E ARE AT THE SAME LOCATION TO SAVE SPACE* 

EQUIVALENCE ( AZ y E ) y ( BZ y F ) 

COMMON /BLKl /TW IRE y TAVE y ALPHA 1 y ALPHA2 y ALPHA3 y BETA y 
1 H 1 y H2 y H3A y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 

COMMON /BLK2/STARTy RANGE 


10 ERFLAG = 0 

MNN =: MN 

THE FOLLOWING STATEMENT IS NEEDED ONLY 
FOR THE EXAMPLE IN THIS REPORT* 

THE MACH NUMBER WAS MEASURED DOWN STREAM OF THE PULSED 
THERMOCOUPLE SITE WHERE THE GAS WAS COOLER* THIS NEXT 
STATEMENT CONVERTS THE MACH NUMBER AT THE LOWER 
TEMPERATURE TO THAT AT THE PULSED THERMOCOUPLE SITE* 

MNN = MN*SQRT(TGAS/MTMP) 

COMPUTE NUSSELT NUMBER* 

^ 15 NU = 1S8*41>K( SORT (MNN*P)K WDIA) )/ 

1 ( ( TGAS**() * 6 ) * ( ( 1 * + * 2*MNN*JK2 ) ** * 25 ) ) 

C 

K1 = WDIA*WDENS*SPHT/(4*)KSIGMA*(El + E2)t:TWF) ) 

C 

K2 = (TGAS*)|c*78))KNU*PSC*3*0()7E-4/ 

1 (WDIA*SIGMA*(E1+E2*TWF) ) 

C 

K3 = K2*TGAS+ ( 1 * -ALPHAG ) * ( TDUCT**4 ) +EGAS* ( TGAS**4 ) 

C 

C 

C COMPUTE ALPHAl yALPHA2yALPHA3yBETAyHl yH2yH3AyH3B* 
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SCALE NUMBERS DOWN BY A FACTOR OF 10**-20 TO 
PREVENT OVERFLOW* 

AZ (K2**2)*(l .0E“-20)*(K2)|oK2)/-4* + 

K3* ( 6A ♦ OE-20 ) * ( K3**2 ) /27 . 0 

BZ = AZ 
Y == 1./3* 

AZ === <SQRT(AZ)>|f(l *OE+J.O) + ( K2**2 ) /2 * ) ** Y 
BZ •” ( K2>K)K2 ) /2 * 0 -• SORT ( BZ ) * ( .1 * GET .1 0 ) 

X == ABS(BZ) 

ERROR IF BZ IS POSITIVE 

IF (X.EQ*BZ) STOP 2 
BZ == X**Y 

EVALUATE ALPHA'S AND BETA* 


Y = AZ-BZ 

ALPHA:!. SQRT(Y)/2* 

BETA == Y + 2*>KK2/SQRT(Y) 

BETA SORT ( BETA )/2* 

X = SO R r ( 2 . K 2 / S R T ( Y •••• Y ) 

ALPHA2 - SORT ( Y ) /2 . f X./.2 * 

A I... P HAS - S a R T ( Y .> / 2 * -• X ./ 2 . 

H:I. = -Kl/( ( ( ALPHA2-ALPHA1 .) >K*2+BETA)(o)f2 .> * ( ALPHA2-ALPHA3 .> ) 
H2 == -K1./( ( (ALPHA3“ALPHA;l. )**2TBETA)!(*2)*(ALPHA3 -ALPHA2) ) 
E = ~ 2 * * ( B E T A * )K 2 ) * < 2 * ♦ A L P H A 1 - A I... P H A 2 - A I... P H A 3 ) 

F 2 * * ( BETA* ( ALPHA 1-ALPHA2 ) * ( ALPHA J. -ALPHAS ) -BETA>K*3 ) 

X = E**2+F**2 
H3A = -K:l,*E/X 
H3B = K;I.*F/X 


TIME == H1.*AL0G(ALPHA2-T) + H2*AL0G < T-ALPHA3 ) + 
H3A*AL0G( (T-ALPHAl )**2+BETA**2) + 
2*0*H3B*ATAN(BETA/(T-ALPHA;l. ) .) + HA 

H4 = CONSTANT TO BE DETERMINED* 

TAVE =••=•■ AVERAGE INITIAL WIRE TEMPERATURE* 

MSTIME IS TIME SCALE FACTOR* 

MSTIME*(TIME INTEGER.) IS TIME SINCE START OF DATA* 

TIME = 0 AT DATA POINT 1=0 

TIME = MSTIME AT DATA POINT 1 = 1. etc* 

SET ERROR FLAG IF ALPHA2 IS LESS THAN TAVE SINCE IT 
WOULD REQUIRE TAKING THE LOG OF A NEGATIVE NUMBER* 


IF (ALPHA2*GT*TAVE) GO TO 40 


no no nooooonnnnnno n noon 


ERFLAG = 6 
RETURN 
C 

40 X = TAUE-ALPHAl 

COMPUTE INTEGRATION CONSTANT H4 BY SETTING T=TA0E 
AT THE STARTING TIME AND SOLOING FOR H4* 

H4 = MSTIME*FLOAT<START) - ( H1*AL0G ( ALPHA2-TA0E ) + 
1 H2*AL0G(TA0E-ALPHA3) + 

1 H3A*AL0G( (TAOE-ALPHAl )*>lc2+BETA**2) + 

;l 2 ♦ 0*H3B*ATAN2 ( BETA » X ) ) 

RETURN 
END 


FUNCTION EOALTM(T) 


PURPOSE 

EOALUATE THE THEORETICAL EQUATION (TEXT EQU* 24) 
FOR TIME AS A FUNCTION OF TEMPERATURE, 

T = INPUT TEMPERATURE (K), 


FUNCTION EOALTM(T) 


REAL TWIRE( 1.000) y TAVE » ALPHAl r ALPHA2 ? ALPHA3 r BETA r 
1 H 1 y H2 r H3A y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 

REAL TyX 

COMMON /BLK 1 /TWIRE y TAOE y ALPHA 1 y ALPHA2 y ALPHA3 y BETA y 
1 HI y H2 y H3A y H3B y H4 y TDUCT y MSTIME y MN y P y MNN y MTMP 


10 X = T-ALPHAl 

EVALTM = H1*AL0G(ALPHA2-T)+ 

1 H2*AL0G(T-ALPHA3)+ 

1 H3A*AL0G( (T~ALPHA1 )**2+BETA>K*2) + 
1 2 . 0*H3B*ATAN2 ( BETA y X ) +H4 

RETURN 
END 
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c 

C SUBROUTINE MATINU 

C 

C PURPOSE : 

C INUERT A 1 OR 2 DEGREE MATRIX. 

C 

SUBROUTINE MATINU ( ARRAY » NORDER ) 

C 

C ARRAY = INPUT MATRIX WHICH IS REPLACED BY ITS INMERSE. 

C NORDER= DEGREE OP MATRIX. 

C 

REAL ARRAY <2*2) f DET y X 
INTEGER NORDERfIfJ 
C 

10 IF (NORDER.EQ. 1 ) GO TO 20 

IF (N0RDER.EQ.2) GO TO 30 

STOP 800 
C 

C CALCULATE INOERSE OF ONE DEGREE MATRIX. 

C 

20 ARRAY(lrl) = l./ARRAYd^l) 

RETURN 

C 

C CALCULATE DETERMINANT FOR SECOND DEGREE MATRIX. 

C 

30 DET = ARRAYd d )*ARRAY(2»2>-ARRAY(1 >-2)*ARRAY(2!> 1 ) 

IF (DET.EQ.O) STOP SOI 
C 

C CALCULATE ADJOINT MATRIX 

C 

X = ARRAY(ItI) 

ARRAY(lfl) == ARRAY (2r2) 

ARRAY (2 1-2) = X 
ARRAY<lf2) == -ARRAY(li-2) 

ARRAY(2d) -ARRAY(2rl) 

C 

C CALCULATE THE INOERSE OF SECOND DEGREE MATRIX. 

C 

DO 50 1=1^2 
DO 40 J=1^2 

ARRAY(ItJ) ••= ARRAYd, J)/DET 
40 CONTINUE 

50 CONTINUE 

C 
C 

RETURN 

END 
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Appendix C 

Gradient -Expansion Method 


This appendix describes the least-squares fit to a 
nonlinear function that uses the gradient-expansion 
algorithm taken from Bevington (ref. 2). The 
objective of the process is to search for the values of 
parameters in the theoretical equation that will 
minimize the sum of the squares of the difference 
between the data points and the theoretical nonlinear 
function. This sum to be minimized is defined as 

m 2 

T,[Yi-Y{Xi)\ (Cl) 

1=1 


Y{X) = YoiX) + 8aj (C4) 

where Yq{X) is the value of Y{X) at the starting point 
for the expansion. Then 


x^= L 

i=l 


Yi-YoiXi) 


j=\ ^ 


(C5) 


where m is the number of data points, T,- is the 
dependent variable, A',- is the independent variable, 
and Y{X) is the theoretical function with unknown 
parameters a^. 

The quantity is regarded as a function of the 
parameters aj of the fitting function Y(X). There are 
m data points {Xj, Yj). The idea is to choose the 
values of the n parameters aj so that is a 
minimum. 

The first approach is to take the gradient of x^ 


Vx^= D 
>=i 


9x^ . 

daj 


(C2) 


where the dj are unit vectors. The gradient of x^ 
gives the direction of the maximum rate of increase 
of x^- We want to increment the parameters from 
some starting value xo so that x^ decreases. Hence 
we write 


8aj = -iVxo)j ^aj 



(C3) 


The Aaj are size constants that must be supplied. The 
parameters aj are incremented by 5a j and the process 
repeated. The minus sign insures that the increments 
are in a direction opposite to the gradient so that they 
are in the direction of most rapid decrease of x^- 
However, the method tends not to work well near the 
actual minimum — it is better further away. 

Another approach is to expand the fitting function 
Y(x) as a first-order Taylor series in the parameters 


We now want to minimize x^ as a function of the 
increments 5aj', so we take d\^/d5ai^ and set it equal 
to zero 


g2[y,-yo(^,)] 


dYoiXj) 

dak 


n m 

= i: 6a, i: 2 

7=1 /=1 


dYf,{Xj) dYQjXj) 
daj dak 


(C6) 


This gives a set of n linear equations for the n 
quantities daj. Define 




1 ^ 

2 dai^ 


m 




Yi~Y(Xi)\ 


dYoiXj) 

dOk 


(C7) 


m 

^jk — 2^ 
1 = 1 


dYo(Yj) dYoiXj) 
daj dOk 


(C8) 


and 


m . 

xl= E |J^i-l^o(^/)| 


(C9) 


thus 


^k = 


n 


y=i 


k=l,2,...,n 


(CIO) 
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This can be put into the form of a matrix equation 
/3 = 6a*a 


where X is an arbitrary parameter that changes the 
method from the Taylor series to the gradient 
method. If X is near zero, the method is the same as 
the Taylor series approach. If X is large, the diagonal 
terms dominate and the equations are essentially 


or 


fi-a ^=5a 


(Cll) 

^j = \bajajj 


where /3 and 6a are column matrices with n elements 
and a is an n-hy-n symmetric square matrix. This 
method tends to work well near the actual minimum 
but poorly far from the minimum. 

By combining the two methods it is possible to 
obtain an algorithm that works well far from the 
minimum and also close to it. To combine the two 
methods, one writes (ref. 7) 


or 

baj = 


\a 


-I^J- 


JJ 


1 

2Xof jj da j 


-1 


2\a 


■JJ 


(C15) 


(3 = a' • 6a (C12) 

where 

<^jk=<^jk for 7^/: (C13) 

and 

ajj=ajj{\+\) forX>:0 (C14) 


which result in the gradient method. 

This technique can be used by starting with an 
arbitrary small value of X, such as 0.001. If the 
computed baj causes to increase instead of 
decrease, the initial guess at the aj is not good 
enough, and is too far from the minimum for the 
second method to work. Then X is increased by a 
factor of 10 and a new set of baj is found. Each time 
X is increased the algorithm is more like just taking 
the gradient, which works well for aj far from 
(oj) . This continues until starts to decrease, at 
whicTi" lime X is divided by 10 at each iteration. By 
this time the minimum will have been found. 
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Appendix D 

Typical Program Input and Results 


This appendix provides an example of data used by 
the computer program. The following data were put 
into the computer program: 

INPUT MACH NUMBER 0.0286 
INPUT PRESSURE IN Pa. 99805. 

INPUT DUCT TEMPERATURE IN DEG. K. 396.0 
INPUT PROBE SHAPE CONSTANT 0.0 
MACH NUMBER TEMPERATURE DEG. K. 415.8 

The following data were put out by the computer 
program: 

GAS TEMPERATURE = 1707.00 K 
FINAL WIRE TEMPERATURE = 1565.79 K 
PROBE SHAPE CONSTANT = 0.850 


The following data were not put out by the computer 
program but may be useful: 

CHISQR = 0.267E + 05 
TAVE = 677.0 
ALPHAl = 1183.2 
ALPHA2 = 1565.8 
ALPHA3= -3932.2 
BETA = 321 8.3 
HI = -0.902 
H2 = 0.259 
H3A = 0.321 
H3B= -0.260 
H4 = 0.072 

The 1000 data points of thermocouple wire 
temperature are shown in the following listing: 
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THERMOCOUPLE TEMPERATURE DATA (K) 


I 

TWIREd) 

1 

674.7 

■:) 

675.7 

3 

676.7 

4 

678.6 

5 

675.7 

6 

677.6 

7 

677.6 

8 

679.5 

9 

679.5 

10 

676.7 

11 

675.7 

12 

675.7 

13 

675.7 

14 

674.7 

15 

677.6 

16 

676.7 

17 

673.8 

18 

675,7 

19 

672.8 

20 

672.8 

21 

670.9 

2'^ 

671.8 

23 

671.8 

24 

671.8 

25 

671.8 

26 

671.8 

27 

672.8 

28 

669.9 

29 

669.9 

30 

671 .8 

31 

672.8 

32 

671.8 

33 

673.8 

34 

673.8 

35 

678.6 

36 

678.6 

37 

678.6 

38 

680.5 

39 

682.4 

40 

680.5 

41 

679.5 

42 

682.4 

43 

681.5 

44 

678.6 

45 

682.4 

46 

683.4 

47 

682 . 4 

48 

681.5 

49 

678.6 

50 

676.7 


I 

TWIREd) 

51 

676.7 

52 

674.7 

53 

676.7 

54 

675.7 

55 

673,8 

56 

674.7 

57 

672.8 

58 

675.7 

59 

671.8 

60 

671,8 

61 

673.8 

62 

673.8 

63 

673.8 

64 

672.8 

65 

671.8 

66 

669,9 

67 

670.9 

68 

674,7 

69 

673.8 

70 

674.7 

71 

675.7 

72 

676.7 

73 

677.6 

74 

677.6 

75 

680.5 

76 

679.5 

77 

680.5 

78 

679,5 

79 

680.5 

80 

680.5 

81 

680.5 

82 

684.3 

83 

684 . 3 

84 

683.4 

85 

683.4 

86 

682.4 

87 

686.3 

88 

685.3 

89 

682.4 

90 

679.5 

91 

678.6 

92 

682.4 

93 

679.5 

94 

681.5 

95 

678.6 

96 

677.6 

97 

678.6 

98 

679.5 

99 

678.6 

100 

678.6 


I 

TWIREd) 

101 

678.6 

102 

678.6 

103 

677.6 

104 

678,6 

105 

678.6 

106 

681.5 

107 

682.4 

108 

684.3 

109 

687.2 

110 

692,0 

111 

694.8 

112 

701.5 

113 

703.4 

114 

708.2 

115 

713.8 

116 

715.7 

117 

717.6 

118 

723.3 

119 

726 . 1 

120 

730.8 

121 

729.8 

122 

734 . 5 

123 

738.2 

124 

738.2 

125 

741.0 

126 

744,8 

127 

746.6 

128 

750.3 

129 

752.2 

130 

754.0 

131 

757 . 8 

132 

762.4 

133 

765.2 

134 

769.8 

135 

773.4 

136 

775 . 3 

137 

779.8 

138 

783.5 

139 

787,2 

140 

789.9 

141 

790.8 

142 

795 . 4 

143 

799.0 

144 

801.7 

145 

807.2 

146 

809,9 

147 

812.6 

148 

814.4 

149 

816.2 

150 

821,6 


I TWIRE(I) 

151 824.3 

152 827.9 

153 828.8 

154 835.1 

155 835.1 

156 837.8 

157 842.2 

158 842.2 

159 844.0 

160 851.1 

161 854.7 

162 856.5 

163 861.8 

164 866.2 

165 867.1 

166 869. 

167 871. 

168 877. 

169 878. 

170 882.1 

171 886.5 

172 889.1 

173 890.8 

174 892.6 

175 896.1 

176 898.7 

177 902.2 

178 903.1 

179 906.5 

180 909.2 

181 912.6 

182 916.1 

183 916.1 

184 918.7 

185 920.4 

186 922.1 

187 925.6 

188 928.2 

189 928.2 

190 931.6 

191 934.2 

192 935.0 

193 935.0 

194 937.6 

195 940.2 

196 944.5 

197 944.5 

198 947.9 

199 950.4 

200 953.0 
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THERMOCOUPLE TEMPERATURE DATA (K> 


I 

TWIRE(I) 

I 

TWIREII) 

I 

TWIRE(I) 

I 

TWIRE( I ) 

201 

957.3 

251 

1070.3 

301 

1176.3 

351 

1256,5 

202 

959.8 

252 

1076.0 

302 

1178.6 

352 

1258.0 

203 

960.7 

253 

1073.5 

303 

1181.8 

353 

1259.5 

204 

964.9 

254 

1077.6 

304 

1181.8 

354 

1259.5 

205 

964.9 

255 

1079.2 

305 

1182.6 

355 

1260.3 

206 

967.4 

256 

1080.8 

306 

1183,3 

356 

1261 . 1 

207 

973.4 

257 

1083.3 

307 

1187.2 

357 

1264.1 

208 

974.2 

258 

1086.5 

308 

1188.0 

358 

1262.6 

209 

975.9 

259 

1088.1 

309 

1189.6 

359 

1264.1 

210 

977.6 

260 

1090.5 

310 

1192.7 

360 

1264.8 

211 

981.8 

261 

1093.8 

311 

1191.9 

361 

1264.8 

212 

981.8 

262 

1093.8 

312 

1193.4 

362 

1267.1 

213 

983.5 

263 

1097.8 

313 

1195.8 

363 

1268.6 

214 

986.0 

264 

1099.4 

314 

1198.9 

364 

1270.2 

215 

986.8 

265 

1103.4 

315 

1199.7 

365 

1270.9 

216 

989.3 

266 

1104.2 

316 

1202.8 

366 

1271.7 

217 

989.3 

267 

1106.6 

317 
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